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Observing a multiple version of the divisor function we introduce a new zeta function which 
' we call a multiple finite Riemann zeta function. We utilize some g-series identity for proving 

the zeta function has an Euler product and then, describe the location of zeros. We study 
further multi- variable and multi-parameter versions of the multiple finite Riemann zeta functions 
and their infinite counterparts in connection with symmetric polynomials and some arithmetic 
quantities called powerful numbers. 
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. 1 Introduction 

° . 

The divisor function dfc(A r ) := J2d\N^ k ^ s a basic multiplicative function and plays an important 

i > i 

role from the beginning of the modern arithmetic study, and in particular, Ufc(iV) appears in the 
Fourier coefficients of the (holomorphic) Eisenstein series Ek+\{r). Usually perhaps, it is not so 
common to regard a^N) as a sort of zeta function. However, in the present paper we treat the 
divisor function cru(N) as a function of a complex variable k = —s. Still, it is clear that when 
a ; oo through LL (or 0, we have U- S {N) «.), the Riemann z eta function, and there 

are at least two interpretations of 

Zh(s) := J>- s = a_ s (iV) 

n\N 

as a zeta function in number theory: 

• Fourier coefficients of real analytic Eisenstein series, 

• Igusa zeta functions. 

Concerning the first, we refer to Bump et al [BCKV] where the so-called "local Riemann hy- 
pothesis" is studied. In case of the real analytic Eisenstein series E(s,t) for the modular group 
SX(2,Z), the iV-th Fourier coefficient is essentially given by 

c n (s,t) := Z l N (2s - l)K s _i(27riVIm(r))e 2 ^ Re M 
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and hence it satisfies the local Riemann hypothesis: if cn{s,t) = then Re(s) = \. On the other 
hand, for the second, the interpretation is coming from the (global) Igusa zeta function £ Igusa (s, R) 
of a ring R defined as 

oo 

C Igusa (s,i?) := Yj # Hom ring ( J R,Z/(m))m- s . 

m=l 

Then, in fact, it is easy to see that 

Z 1 N (s) = ( l ^(s,Z/(N)). 
Thus, the purpose of the present paper is initially to study the function defined by the series 
(1.1) Zn{s):= (^2---n m )- s . 

ni\n2\---\n m \N 

We call Z™(s) the multiple finite Riemann zeta function of type N. We then study several basic 
properties of Z™(s) such as an Euler product, a functional equation and an analogue of the Riemann 
hypothesis in an elementary way by the help of some g-series identity. We also study the limit case 
Z Z( S ) '■= limAf-KXD Z%(s)(= Zjp(s)) = J2 ni \n 2 \-\n m (. n i n 2 ■ ■ ■ n m )~ & '■ Here and throughout the paper 
we consider the limit N — > oo as the one through factorials. 

Moreover, we generalize this zeta function in two directions. The first one is to increase the 
number of variables. We prove that the Euler product of a multi-variable version of the zeta function 
is expressed in terms of the complete symmetric polynomials with a remarkable specialization of 
variables (Theorem 3.2). The second one is further to add parameters indexed by a set of positive 
integers. For general parameters it seems difficult to calculate an explicit expression of the Euler 
product using symmetric functions with some meaningful specialization of variables. However, 
when we restrict ourselves to the one variable case, under a special but non-trivial specialization 
of parameters, we show that the corresponding multiple zeta functions are written as a product of 
Riemann's zeta functions ((cs) with several constants c's determined by the given parameters and 
the Dirichlet series associated with generalized powerful numbers (see Section 4 for the definition). 
Moreover, we determine the condition whether the Dirichlet series associated with such generalized 
powerful numbers can be extended as a meromorphic function to the entire complex plane C or not 
(see Theorem 4.8 and its corollary). As a consequence, the most of such Dirichlet series are shown 
to have the imaginary axis as a natural boundary. The result is a generalization of the one in [IS]. 

In the final position of the paper, we make a remark on the relation between the multiple finite 
Riemann zeta functions and the number of isomorphism classes of abelian groups. In particular, 
using the Tauberian theorem, for a given positive integer m, we prove that the asymptotic average 
of the isomorphism classes of abelian groups of order n which are given by the direct sum of p- 
groups A p such that p m A p = 0, equals £(2)£(3) • • • ((m) (when n tends to oo). We make also a small 
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discussion on an analogous notion of (holomorphic) multiple Eisenstein series which are defined via 
these multiple finite Riemann zeta functions as coefficients of its Fourier series (or related to the 
so-called Lambert type series). 

Throughout the paper, we denote the sets of all integers, positive integers, non-negative integers, 
real numbers and complex numbers by Z, Z>o, Z>o, R and C respectively. 

2 Multiple finite Riemann zeta function 

In this section we prove the fundamental properties of the multiple finite Riemann zeta function 
Z™(s) defined by (1.1) and make a discussion on some related Dirichlet series. Namely we first 
show the following theorem. 

Theorem 2.1. Let N be a positive integer. 

1. Euler product : 

„ ™ 1 _ v s(ord p N+k) 

z n{s)= n n i_ p -.* ■ 

p : prime k=l 

where oid p N denotes the order of p- factor in the prime decomposition of N. 

2. Functional equation : 

(2.2) Z%(-s) = N ms Z%(s). 

3. Analogue of the Riemann hypothesis : All zeros of Z™(s) lie on the imaginary axis Res = 0. 
More precisely, the zeros of Z^{s) are of the form s = ( or( j w+ky\ogp f or ^ = 1> • • • > m > 
and n G Z \ {0}. Consequently, the order Mult m (n,p, k) of the zero at s = ( or d Ar+fc)iogp ^ s 
given by Mult m (n,p,k) = #{(/,j) ; 1 < I < m, j G Z\{0}, (ord p N + k)j = (oid p N + l)n}. 

4- Special value : When n is a positive integer one has Z^(—n) £ Z. 

For the proof of the theorem the following lemma is crucial. 

Lemma 2.2. Let m be a positive integer. Then 

1. For any integers I > it holds that 

( 2 -3) E 

where is the q-binomial coefficient defined by [^] = I^jLiU — q n+l ~i) / '(1 — qi). 



m — 1 + d 




m + 1 


m — 1 






q 


m 
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2. For \x\ < 1, \q\ < 1 it holds that 

CO 

(2-4) E 



(2=0 



m + d 
m 



x 



m 1 



k=0 



q k x 



Proof. We prove the formula (2.3) by induction on I. When I = the formula (2.3) clearly holds. 
Suppose that it holds for /. Then we see that 



l+i 

E 



m — 1 + d 




m + I 


+ 


m + I 


m — 1 






m — 1 


i 


m 


Q 



J+l 



m— 1 



ii( i - fr 1 n a - ^ m+fc ) {(i - + (i - 9 m )^ +i } 

k=l k=l 
m m 

na-^IK 1 -^) 



k=i 



k=i 



m + l + l 

m 



whence (2.3) is also true for I + 1. 

The second formula (2.4) can be proved in the same manner by induction, but on m. (It is also 
obtained from the so-called (/-binomial theorem. See, e.g. [AAR]). This completes the proof of the 
lemma. □ 

Proof of Theorem 2.1. The functional equation (2.2) is easily seen from the definition. Actually, 
we have 

E (nm 2 ...n m y = N™ £ (^■■■lf) S = NmSZ N^)- 
ni|ri2|---|n m |A r rci|n 2 |---|n m |iV 

Next we prove the Euler product (2.1). We first show that Z^}(s) is multiplicative with respect 
to N. In fact, if we suppose N and M are co-prime, then we observe 



^2 ( n i n 2 • • • n m ) s 

n 1 \n 2 \---\n m \NM 

E ( c i c 2 ■ ■ ■ c m )~ s E {did 2 ---d ri 



Z™(s)Zm(s), 



Cl\c 2 \-\Cm\N 



di|ci 2 |-|rf m |M 



because every divisor n of NM is uniquely written as n = cd where c\N and d\M. By means of 
this fact, in order to get the Euler product expression (2.1) of Z^(s) it suffices to calculate the case 
where TV is a power of prime p. In this case, one proves the expression 



(2.5) 



rn _ { l +k)s 

z%{s) = n — 



k=l 



P 



-sk 
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by induction on m as follows: It is clear that (2.5) holds for m = 1. Now, suppose (2.5) is true for 
m — 1. Then we see that 



P 

0<jl<-<jm<l 



m_1 1 _ p -(d+k)s 



m — \ + d 
d 



p~ ds . 



= E ^Mp-*- = E n i-^ p ~ ds = s 

Therefore the assertion follows immediately from the formula (2.3) in the lemma above. This proves 
(2.5), whence the Euler product for Z™(s) follows. 

Using the Euler product (2.1), we observe that the meromorphic function Z™{s) may have zeros 



at 



each s = ( or( j iv+fc) log p for fe = 1, . . . , m, p\N and n G Z. Note however, since 



(2.6) ^ ( o)= n r p ; >°< 

p : prime 

s = is not a zero of Z™(s). Suppose now 

2mn 2nim 



(ordp N + k) log p (ordg TV + I) log g 

holds for some 1 < k, I < m, p, q\N and n, j G Z \ {0}. Then it becomes p( ord f jV+fe ^' = g( ord " JV+ ^ n . 
This immediately shows that p = q and (ord p N + k) j = (ord p N + Z)n. Hence the order of zero at 
s = (ord W+k)Togp 1S ac t uan y given by Mult m (n,p, k). Obviously, one has Mult 1 (n,p, 1) = 1. 

The last claim about the special value of Z™(s) is clear from the definition. This completes the 
proof of the theorem. □ 

Remark 2.1. Note that Z™(s) = Z l pm (s) from (2.5). □ 
Letting TV — > oo (or N — > 0) in Theorem 2.1, we have the following 

Corollary 2.3. Define a function Z™(s) by 

Z™{s):= E (n 1 n 2 ---n m )- a (Res > 1). 

m\n2\---\n m 

Then we have Z™(s) = IJfcLi C(ks). Here ((s) is the Riemann zeta function. □ 

We next consider a Dirichlet series defined via the multiple finite Riemann zeta functions. 
Corollary 2.4. Retain the notation in the corollary above. Define a Dirichlet series C m ( s ) by 

oo 

rw :=E z « m ( s ) re ~ s - 

71=1 
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Then we have 

m+l 

(2.7) ( m (s)=Hask) = ZZ + \s). 

k=l 

Proof. Since Z™(s) is multiplicative relative to n, the Dirichlet series ( m (s) has an Euler product 
of the form 

C m (s)= H fczytfp-A. 

p : prime \ (=0 / 

Hence the first equality follows from Lemma 2.2. In fact, 

00 oo m -s(l+k) m + 1 -I 

E ^Mp- = E II fear-p- 1 - II rr^i 

Z=0 1=0 k=l k=l 

when Re s > 0. The second equality is obvious from Corollary 2.3. This proves the corollary. □ 

Remark 2.2. Let g(n) be the number of the isomorphism classes of the abelian groups of order n. 
Then we have (see, e.g. [Z, p. 16]) 

00 00 
lim ( m (s) = fl C(ks) = ^g{n)n- s (Res > 1). 

k=l n=0 

Observing the analytic property of Z'^ l (s), we have much precise information about the isomorphism 
classes and their asymptotic averages. For details, see Section 5. □ 

We now slightly generalize the discussion above. Put 



in "*. 



CM) :=J2 Z n(s)i 
n=l 

It is clear that C m ( s ) = C m ( s ) s )- Then, by the same reasoning, ( m (s,t) has the Euler product 

(00 \ m 

1=0 / k=0 

Actually, it follows from the formula (2.4) that 



00 / m 



1 _ p -s(k+l) \ rn 1 



Z^j \ II I _ v -sk P ) II 1 , - i 

1=0 \k=l F J k=0 F 

for Res > and Ret > 0. This is a generalization of the formula (2.7). One may consider C m (s,t) 
as a zeta function for two variables. We give here a few examples: 

Example 2.1. The following is quite well-known. For Ret > k + 1, we have 

00 00 
CH-k, t) = J2 Z'ni-kh-t = c-k^n-t = C(t)C(t - k). 

n=l n=l 

□ 
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□ 



Example 2.2. From (2.6), for Ret > 1, it follows that 

cw) = fxwn-* = f; n ( m+ " dpn )n' t = at) m+1 - 

n=l n=lp: prime 

Example 2.3. Let I be a positive integer. For Ret > 1 — I, we have 

771+ 1 

c m (i,t+i) = i[ at + ik). 

k=l 

As m -> oo, linim^oo ( m (l, t + /) = flfeLi C(* + ^) is the higher Riemann zeta function studied 
in [KMW] (see also [KW]). It should be noted that the higher Riemann zeta function Czoo(s) := 
rife^=i C( s + ^k) possesses a functional equation. □ 

3 Multivariable version 

We generalize the definition (1.1) of the multiple finite Riemann zeta functions Z™(s). For 7 = 
(71, ... , 7 m ) G Z> and N € Z >0 , we define 



(3.1) Z^(t 1; ...,t m ) := Y, ' 

nT\-\n^\N 



This is multiplicative with respect to N. We also notice that Z^(s) = Zj l (s, . . . , s) when 7 

rrt 

(l m ) = (1, . . . , 1). We can prove the following lemma in a similar way to Theorem 2.1. 
Lemma 3.1. For I > 0, define a function Gj(qi, q%, . . . , q m ) by 



(3.2) G]( gi ,...,q m )= Y, Vi 1 ---Vt, 

«>Ai>A 2 >->A m >0 

where the sum is taken over all partitions A = (Ai, A2, . . .) of length < m such that Ai < I and 
7j|Aj for 1 < j < m. Then the Euler product of Zjj(ti, . . . , t m ) is given as 

zj f (t 1 ,...,t m ) = J] g^nCp-'S-.p"*-)- 

p : prime 

□ 

Remark 3.1. Let /(ni, . . . , n m ) be a function defined on G Z™ which is multiplicative with respect 
to each variable rij. Then, in general, we can define a multiple zeta function by 

z nU)= Y f(ni,...,n m ). 
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Actually, one can show that Z 1 N {f) is multiplicative with respect to N, whence has the Euler 
product as (/) = \[ p . primc Zj oldp N (f). □ 

We first look at the simplest case where 7 = (l m ). We abbreviate respectively Z N '(ti, . . . , t m ) 
and \q 1 , . . . ,q m ) to Z%{ti, . . . ,t m ) and Gf{qi, . . .,q m ). 

Theorem 3.2. Let hj(x\, . . . ,x m ) be the j-th complete symmetric polynomial defined by 

hj (x\ 1 ■ ■ ■ , Xfn) = ^ ^ X\ X^ . 
iiH Hm=j 

T/ien we have 

1 

(3.3) . . . ,q m ) = ^h j {q 1 ,q l q 2 , . . . ,q x q 2 - ■ ■ q m ). 

3=0 

In particular, we have 

m 1 

(3.4) GZ(qi, ...,q m ):= Km G^qi, ,q m ) = U z ■ 

l-^oo ^ 1 - qiq 2 ■■■qk 

Proof. Since 



U ( „ „ „ „ „ „ \ _ \ " il+l2 J i Mm „«2H Mm S m 

nj [qi, qiQ2, qiq2- ■■ q m ) — 2_> q i q ? •" < ? m 

•il+«2H Mm=j 

E ?W > 

j=Ai>A 2 >->A m >0 

the first formula (3.3) is clear from the definition (3.2). The second formula (3.4) follows from (3.3) 
together with the fact that the generating function of complete symmetric polynomials is given by 
T,jLohj(xi, . . . ,x m )z J = IlfcLiU - Khz)" 1 (see [M]). □ 

As a corollary of the theorem, we obtain the Euler product of Z™(t±, ... ,t m ). 

Corollary 3.3. Let m and N be positive integers. Then we have 

Z%( tl ,...,t m )= £ n^-n-*™ 

n m |-"|m|JV 

(ordp N 
j=0 

Further, for Retj > 1 (1 < j < m), it follows that 

m 

ZZ(h, . . . ,t m ) := E ^ h ---n-^ = l[C(t 1 +t 2 + ---+t k ). 

n m \---\n 1 k=l 



□ 
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The following lemma gives a recurrence equation among G](qi,..., <7 m )'s. 
Lemma 3.4. For 7 = (71, ... , j m ) £ Z™ 0; we have 

n=0 

.Here |_a^J is i/ie largest integer not exceeding x. 
Proof. Observe that 

LV71J 

G7(5i,...,ftn)= E ? 71Tt E 9 2 A2 ---^= E ^rGfe-'^fe,...,^)- 

0<n<i/7l 7in>A 2 >--->A m >0, n=0 

7j|Aj (2<i<m) 



□ 



The following lemma shows that it is enough to study the case where 71 , . . . , 7 m are relatively 
prime. 

Lemma 3.5. For 7 = (dc\, dc 2l . . . dc m ), we have 

G](q 1 ,...,q m ) = G^ c -\qf,...,q d m ). 

Proof. It is straightforward. □ 

Let us calculate several examples of G' y (qi, . . . ,q m ) for special parameters 7 which are rather 
non-trivial. 

Example 3.1. Let 7 = (c, 1) and calculate G < f ,1 \q\,q2)- Putting d = \ l/c\, we see by Lemma 3.4 
that 

G\ c > 1 \ qi , q2 ) = ^ q <TGU q2 ) 

n=0 

d 1 n Cn + 1 1 / 1 r , C ( d+1 ) 1 (r, r, \c(d+l) \ 

(3.5) = Vaf n ^^ = ^—(1^1 i-(gig 2 )i 

I-Q2 1-aA l-a? 92 l-(q,q 2 Y 



?1=0 



(1 - qf - (g lg2 ) c ) - g2 (l - gf)(l - (gig 2 ) c(d+1) ) 
(l- g2 )(l-^)(l-( gi g 2 ) c ) 

If we take a limit / — > 00, we obtain 

1 - Q2 + ^2 - <7i<? 2 



Coo Wl,92j (1 _ ?2)(1 _ g c )(1 _^ 2C) 
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Therefore we have 



z£ 1 \h,t 2 )= £ W 2 

712 |nf 

= C(i 2 )C(rfi)C(c(ii + t 2 )) H (1 - p-* 2 + p- ct1 -' 2 - p-*^] 



p : prime 



for Retj > 1 (j = 1,2). Thus we may have various possibility of Euler products of the form 
U P : prime ff(p- a ,P - *), where T) G 1 + S-C[S,T] +T-C[S,T] is arising from Z& 1) (t 1 ,t 2 ). □ 

Example 3.2. We calculate G OQ d,c ' 1 \q, q, q). If we set q± = q 2 = q in (3.5), then we have 

= (i-g c(d+1) )(i - g + g c - g cd+c+1 ) 



(l-g)(l-g2c) 
It follows from Lemma 3.4 that 



n=0 n=0 V y ^ y ' 

1 f 1 - q + g c <f (1 + q c ) q 2c+1 



(1 - - g 2c ) 1 1 - q cd 1 - q 2cd 1 - 9 3cd 

In particular, if we put d = 1 then 

1 - q + q c - q c+1 + q 2c+1 - q 3c + q 3c+1 - q 4c 1 - q + q c - q c+1 + q 2c 



Gfe c>1) (9,9,9) 



(1 - g)(l - g 2c ) 2 (l - g 3c ) " (1 - g)(l - g 2c )(l - g 3c ) ' 

Hence we obtain 

Z^ c ' 1 \s)= Y, (nln^y* = asKilcsnScs) [J (1 - p~ s +p-» - p-< c+1 ) s + p- 2cs ) 

n 3 |n^|nf p : prime 



□ 



4 Multiple zeta functions and powerful numbers 

In this section, we study Z2o{s) for 7 = (k, k, . . . , k, 1) in connection with a certain generalized 
notion of powerful numbers. 

Let us first recall the definition of powerful numbers. A positive number n is called a /c-powerful 
number if ord p n > k for any prime number p unless ord p n = (see, e.g. [IS], [I]). Extending this 
we arrive at a new notion, l-step k-powerful numbers; a positive integer n is said to be an /-step 
/c-powerful number if n satisfies the condition that ord p n = 0, k, 2k, ...,(/ — l)k or ord p n > Ik 
for any prime number p. Clearly, if n is an /-step /c-powerful number, then n is again a .7 -step 
/c-powerful number for each j (1 < j < I). In particular, l-step /c-powerful numbers are nothing 
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but the usual /c-powerful numbers. Note also that every natural number is an /-step 1-powerful 
number for any /; this agrees with the claim for k = 1 in Theorem 4.8 below. 

As an example of /-step fc-powerful numbers, for instance, the first few of 2-step 2-powerful 
numbers are 1, 4, 9, 16, 25, 32, 36, 49, 64, 81, 100, 121, 128, 144, 169, 196, 225, 243, .... We note that in 
general, an /-step fc-powerful number n has the following canonical representation: 

n = a k ia f ■ ■ ■ af~ l)k x m, 

where a±,- ■ ■ ,ai are square-free, m is /A;-powerful and these satisfy gcd(ai, • • • , a;,m) = 1. Note 
that here a fc-powerful number m is uniquely expressed as m = bib% +l ■ ■ ■ b 2 ^~ l if we stipulate that 
^2, • • • bk are all square- free. 
Let us put 



fk,l( n ) :-- 



1 if n is an /-step fc-powerful number, 
otherwise 



for a positive integer n. We define also Fkj(s) := YHn=i fk,i( n ) n ~ s - This arithmetic function fk,i{n) 
is multiplicative with respect to n. Note that Fij(s) = ((s) for any /. We show that Z^o(s) is 
represented by the product of the Riemann zeta functions times Fkj(s). Namely, we have the 

l 

Theorem 4.1. Let k, I be positive integers and put 7 = (k, k, . . . , k, 1). Then we have 

1+1 

(4.1) Zl(s)= Y, (4---™Nmr = ^(s)n^) (Re S >l). 

ni +1 \nf\-\n'l j=2 

Remark 4.1. When k = 1, this theorem gives Corollary 2.3. □ 
To prove the theorem, we prepare the following two lemmas. 



Lemma 4.2. Let 7 = (k, k, . . . , k, 1). Then we have 

'+ 1 l +l 1 1 _ „ + Jk+i _ Ml+i) 
GUI) = GUtt^Ta-) = II T—fk ■ " \_ q " 

i=i 
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Proof. By definition, it follows that 

oo 

G (fc,...,fc,l)( g ) = ^ gAi+-+Ai+A i+ i = y^ q n q M+-+M 

Ai>->A ; >A i+1 >0 n=0 Ai>->A ; >n 

k\\j (l<j<l) klXj (l<j<l) 

oo fc 

_ ^ q Xi+-+\i _|_ ^2 ^2 q ak+b ^2 q x i+-+ x i 

Ai>--->A ; >0 a=0 6=1 A!>--->A ; >afc+6 

(l<i<0 fe|Aj (l<j<i) 

oo k 

= ^2 g fc(m+ " ,+w) + q ak+b y q k{m+-+m) 

/ii>--->^t ; >0 a=0 6=1 Mi>--->M;>a+l 

oo k 

= G l OQ (q k ) + ^2q ak ^2q b ^ q k{(vi+a+i)+-+(v l +a+i)} 
a=0 6=1 ^i>--->j/ ; >0 

= gUq") w fe (fy (,+1)o> ) (jy) ( E 

Va=0 / \6=1 / \z/i>-->^>0 



Since G^g) = nj=i(l - Q j )~\ we have 

G ( k ,...,k,i ){ ) = tj _J_ ( 1 + g*' 

{q) 1 = 1 1 - gi* V 1 - g fc (' +1 ) 1 - « / 



1 (l-g)(l-^( /+1 ))+^ fe+1 (l-^) 



II T gJ* i - g 

1 1 _ q _|_ g'fc+l _ 



l+l 

This proves the assertion. □ 
Lemma 4.3. We have 

g + \_ = (l-q k )[l+q k + q 2k + --- + q^ k + q lk y j < 

Proof. The calculation is straightforward. Actually, we have 

(oo 
i + q k + g 2fc + • • • + q {l ' l)k + q lk E ? 

= r ^{(l-9 ,fc )(l-5) + 9 ,fc (l-9*)} x _ g 

□ 



1 _ g _)_ gifc+i _ qk{l+i) 
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Proof of Theorem 4-1- It follows from Lemma 4.2 and Lemma 4.3 that 

z&~> k > i \s)= n Gt- Ai \p- s ) 



p : prime 

oo 



- 11 11, ' 

p : prime \ j=2 

(i + p~ ks + p - 2ks -i i- p-Q- 1 )^ + p~ iks + p-( ik + i ) s _i ) 



X 

p : prime 
l+l 

= UC(jks)-F k ^s). 

3=2 

This completes the proof of the theorem. □ 

Now we determine the condition if the Dirichlet series can be meromorphically extended 

to the whole complex plane C. We recall the following Estermann's result [E] (see [K] for a 
generalization). A polynomial f(T) 6 1 + T ■ C[T] is said to be unitary if and only if there is a 
unitary matrix M such that f(T) = det(l - MT). 

Lemma 4.4. For a polynomial f(T) G 1 + T ■ C[T], put L(s, f) = \[ p . primc f{p~ s ). Then 

1. f(T) is unitary if and only if L(s,f) can be extended as a meromorphic function on C. 

2. f(T) is not unitary if and only if L(s,f) can be extended as a meromorphic function in 
Re s > with the natural boundary Re s = ; each point on Res = is a limit-point of poles 
of L(s,f) in Res > 0. 

□ 

Since 



oo 



Fk,i(*)= n ii+p- ks +p- 2ks +---+p- (i - i)ks +p- iks Y,p~ 

p: prime \ j=0 

= as)aks) n (i-p-*+p-(»+i)»-p-*c+i)«) 

p : prime 

by Lemma 4.3, we have only to see whether the polynomial Gk t i{T) := 1 — T + T lk+1 — T k ( l+1 ^ 
is unitary or not. The polynomial Gt,i{T) can be expressed as G k ,i(T) = (1 - T k )H K i{T) with 
H Kl {T) := 1 + (T k - T) J2 l f=o Tkj - 

Proposition 4.5. Let k and I be positive integers. The polynomial G^(T) is unitary if and only 
if k = 1,2. 
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In order to prove this proposition, we need the following two lemmas. 

Lemma 4.6. Let k be a positive integer such that k > 3. Then the unitary root a (i.e. \a\ = 1) of 
the polynomial G k ^(T) must satisfy a k = 1 or a k ~ 2 = 1. 

Proof. Let a = e 2ni6 / 1 (0 e R) be a unitary root of G k ,i(T). Since G M (T)/(1-T) = l+T ifc+1 (l- 
T fc - 1 )/(1-T), wehavel+a /fe+1 (l-a fc_1 )/(l-a) = so that |(1 - a*" 1 )/^ - a)| = |a^ fc+1 )| = 1. 
Hence we see that Re q^" 1 = Re a, that is, cos 2ir(k — 1)0 — cos 2n6 = —2 sin irk9 sin n(k — 2)6 = 0. 
Thus we conclude that either k6 € Z or (A; — 2)6* G Z. This proves the lemma. □ 

Lemma 4.7. Zei k be a positive integer such that k > 3. Suppose that a complex number a satisfies 
a k - 2 = 1. T/ien G'^(a) / 0. 

Proof. Since G£,(T) = (Ik + l)Z/cT^ fc - 2 )+ 2Z - 1 - (ifei + k)(kl + k- l)T( fc - 2 )( ,+1 ) +2, J if we assume that 
a satisfies G£,(a) = and a fe ~ 2 = 1, we have G'^(a) = (/fc + l)/^- 1 - (fc/ + fc)(fc/ + A;-l)a 2 ' = 0. 
This immediately shows that a = (^[jp^pfcirq i which contradicts a fc ~ 2 = 1. Hence the assertion 
follows. □ 

Proof of Proposition 4-5. Let I be a positive integer. Since the unitarity of G k j(T) and that of 
Hk,i(T) are equivalent, it suffices to check the unitarity of H k j(T). It is clear that H k ,i(T) is a 
unitary polynomial when k = 1, 2. Actually we have -Hi,z(T) = 1 and i? 2 ,/( r ) = 1 ~ T + t2 ' +1 ~ 
T 2Z+2 = (i _ + T 2l+1 ), which are indeed unitary. 

Suppose k > 3. If Gk,i{T) is unitary, then every root of H k j(T) satisfies a fc = 1 or a fc ~ 2 = 1 by 
Lemma 4.6. However, if a k = 1 we immediately see that H k j(a) = 1 + (1 — a)l. Thus, H k j(a) can 
not be because of the unitarity of a. Thus any root of H k ^(T) must satisfy a k ~ 2 = 1 and a fc / 1. 
By Lemma 4.7, the multiplicity of these roots of H k j(T) is at most 2. Since H k j(T) is assumed to 
be unitary and H k j(\) / 0, it follows that 2(k — 3) > deg H k i(T) = Zfc. This is possible only when 
/ = 1. 

Therefore it is enough to prove that H ki \(T) is not unitary for k > 3. We put H k (T) = 
H kA (T) = 1 - T + T k for simplicity. 

Suppose that k is an odd integer such that k > 3. Then H k (T) has a real root in the interval 
(-2, -1) since H k (-1) = 1 > and # fc (-2) = 3 - 2 fc < 0. This implies the polynomial H k (T) is 
not unitary. 

Thus, only we have to consider is the case where k is even and k > 4. Suppose that H k (T) 
is unitary and let e* 61 (— 7r < 6 < it) be its unitary root. Then we see that 6 satisfies the 
equations cos k9 = cos 6 — 1 and sinkO = sin 0. Since 1 = sin 2 kO + cos 2 k9 = 2 — 2cos0, we 
have 9 = ±7r/3. Further, since 1 = cos 2 6 + sin 2 6* = (cos kO + l) 2 + sin 2 k9 = 2cos k6 + 2, we 
have cos(fc-7r/3) = —1/2. Hence we see that either k = 2 or 4 (mod 6). On the other hand, since 
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1 = (cos 9- cos k0) 2 + (sin 9- sin k9) 2 = 2 - 2 cos(/c - 1)9, we have cos((fc- 1)tt/3) = 1/2. It follows 
that either k = or 2 (mod 6) holds. Thus we have A; = 2 (mod 6). 

Now we show that every unitary root of H^iT) is simple. If we assume that (3 is a multiple root of 
H k (T), it follows that (3 k -f5+ 1 = and kj3 k ~ l -1 = 0. Then we have |/3| = fc-V(*-i) and /? = 
by the second equation. On the other hand, by the first equation, we obtain 1 = (3 — fi k = (k — l)/3 k 
so that \/3\ = (k- l)- 1 / fc . Therefore we have k k = |/3p fe(fc_1) = (k - l)*" 1 , but this contradicts to 
the assumption of the unitarity of H).(T). This completes the proof of the proposition. □ 

Finally we obtain the following theorem, which gives actually a generalization of the result in 
[IS] concerning the powerful numbers. 

Theorem 4.8. Let k and I be positive integers. When k = 1,2 we have 

„ ( , n , P ( , C(2s)(((2l + l)s) 
Fi,iW=CW, F2,i(s)= C[2{2l + 1)s) ■ 

When k > 3, ^/(s) can be extended as a meromorphic function in Res > and has a natural 
boundary Re s = 0. 

Proof. The assertion is clear from Lemma 4.4 and Proposition 4.5. □ 
As a consequence we have the following result. 

(k k 1) / 

Corollary 4.9. Zet Z£ ,m "' * ; (s) 6e as in Theorem 4-1- Then we have for k = 1, 2 

C((2Z + l)s) ^ 
C (2(2/ + l)s)ll' 



and for k > 3 ^/ie function z!£''"' k ' 1 \s) can be meromorphically extended to the half plane Res > 
with a natural boundary Res = 0. □ 

5 Closing remarks 

We make two small remarks; the first one is to describe the relation between a multiple finite 
Riemann zeta function and the theory of elementary divisors, and the second one is to examine a 
multiple Eisenstein series defined via the multiple finite Riemann zeta function. 

• The isomorphism classes of abelian groups A of order n are parametrized by the map A from the 
set of all prime numbers to that of partitions such that n = FJ p . primc p\ x( - p ^ and 

e(\(p)) 

A - %/p Xjip) z, 

p : primc j=l 
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where |A(p)| and £(X(p)) are the size and the length of the partition \(p) = {\j{p))j>i respectively. 
The multiple finite Riemann zeta function Z^(s) is expressed also as Z™(s) = Yl n \N m 9N( n ) n ~ s • 
Here g™(n) denotes the number of isomorphism classes of abelian groups of order n, parametrized 
by A such that X\(p) < m and £(X(p)) < ord p N for all p. It is clear that g^}{n) is multiplicative 
with respect to n and N. If we put g^{n) := limjv->oo g™(n), then g^(n) equals the number of the 
isomorphism classes of abelian groups A of order n which is the direct sum of p-groups A p such 
that p m A p = for p\n. 

We now just mention about the asymptotic average for g^{n) and Z™(o~) (a G R) with respect 
to n. Thus we need the Tauberian theorem below (see, e.g. [MM]). 

Lemma 5.1. Let F(t) = YlnLi ^n^ 1 be a Dirichlet series with non-negative real coefficients which 
converges absolutely for Re(i) > 0. Suppose that F(t) has a meromorphic continuation to the region 
Re(t) > (3 with a pole of order a + 1 at t = for some a > 0. Put 

c:=— t lim(i- 0) a+1 F(t). 

Then we have 

J> n = (c + (l))^(logx) Q 

n<x 

as x — ► oo. □ 
Using the lemma above, we have easily the following facts: Let m be a positive integer. 
1. 



(5-1) ^C(n) = (C(2)C(3)---C(m) + o(l)) 5 



n<x 



as x — > oo. In other words, the asymptotic average of g^(n) with respect to n is given by 
C(2)C(3)---C(m). 

2. For a fixed a > 0, 

(5.2) Yl = (C(* + l)C(2a + 1) • • • C(w + 1) + o(l))x, 

(5.3) Yl Z ni-°) = (C(* + 1)C(2<7 + 1) • • • C(m<7 + 1) + o{l))x l+m % 

n<x 

(5.4) yj zno) = e n ( ordp : + '") = ih. + "«) *< 1o ^>'" 



n<x n<x p : prime 

as x — ► oo. 
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Actually, since ((s) has a single pole at s = 1 and Res s= i ((s) = 1, it follows the first formula 
(5.1) from (2.7) and Lemma 5.1. Next, fixuei By (2.8), we have C m (a,t) = Y^=i Z n^)n~ l = 
llfclo C(t + ka). This shows that the abscissa of absolute convergence of the Dirichlet series C m (°"; *) 
is given by t = max{l, 1 — ma}. Hence the remaining formulas follow similarly. 

Remark 5.1. Put g(n) := \im m ^ ^ g^(n). It is well-known that J2 n <x di n ) = ( T\T=2 ((k)^J x + 

0(y/x). □ 

Remark 5.2. Since Z^(0) = d(n) := Yld\n ^ follows from (5.4) that ]Cn<x <^( n ) ~ xlogx. It is 
also well-known [Z] that there exists a constant C such that X}n<:r = x 1°S x + Cx + 0{yfx) in 
an elementary way. See, e.g. [A]. □ 

• We define a multiple Eisenstein series with parameter s of type m by 

oo 

ET(q) = Y, Z n(l-s)q n . 

n=l 

We sometimes write E™(t) instead of E™(q) when q = e 2ntT with r G C, Im(r) > 0. It is obvious 
that E\{q) is (essentially) the usual holomorphic Eisenstein series of weight k. In this remark we 
make an experimental study of this multiple Eisenstein series. Here we assume that m = 2; we 
shall make a detailed study of these series in somewhere else. 
First we observe the following simple relation. 

Lemma 5.2. We have 

oo oo 
1=1 N=l 

Proof. The calculation is straightforward. Actually, one has 



oo oo oo oo oo 



Z=l AT=1 Z=l m=l ra=l /=! j=l n \j 



oo oo 



= EEE^V" = E ^(-^ = ^+^)- 

This proves the assertion. □ 

Recall now the Fourier expansion of the holomorphic Eisenstein series £ , fc+i(r) of weight k + 1 
with /c being odd; 



E k+ i(r) = \ £ (cr + d)-*" 1 = 1 + £ J>r + n) 

(c,d)=l V ; m=lneZ 



-it 



1 (27Tz) fc+1 A , . n ( , 1 (2vri) fc+1 . . ; 

1 + C (fc + i) « g> (n)g ( =1 + C ( fc + i) fci 
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Taking fc-times derivative of E^ +1 {t), we have easily 

,v ' n=l 

Hence by Lemma 5.2, we immediately obtain the expression of £'| +1 (r) when k is odd: 

Remark 5.3. There is another expression of £?| +1 (t) for k being odd as follows. 

(2k)\ 



V ; /=1 (c,d)=l 

OO 



□ 



• It is clearly interesting to make an extensive study of the multi-variable and multi-parameter 
version of multiple finite Riemann zeta functions from the symmetric functions point of view. 
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